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ABSTRACT 

Context. Silicate carbon stars are characterized by oxygen-rich circumstellar environments as revealed by prominent silicate emis- 
sion, despite their carbon-rich photospheres. While the presence of a circumbinary disk or a disk around an unseen, low-luminosity 
companion is suggested to explain the peculiar dust chemistry, the origin of silicate carbon stars is still a puzzle to date. 
Aims. We present multi-epoch high-angular resolution observations of 22 GHz H2O masers toward the silicate carbon star EU And to 
probe the spatio-kinematic distribution of oxygen-rich material. 

Methods. EU And was observed at three epochs (maximum time interval of 14 months) with the Very Long Baseline Array (VLBA). 
Results. Our VLBA observations of the 22 GHz H2O masers have revealed that the maser spots are distributed along a straight line 
across ~20 mas, with a slight hint of an S-shaped structure. The observed spectra show three prominent velocity components at Vlsr 
= -42, -38, and -34 km s"', with the masers in SW redshifted and those in NE blueshifted. The maser spots located in the middle 
of the overall distribution correspond to the component at Vlsr = -38 km s"', which approximately coincides with the systemic 
velocity. These observations can be interpreted as either an emerging helical jet or a disk viewed almost edge-on (a circumbinary or 
circum-companion disk). However, the outward motion measured in the VLBA images taken 14 months apart is much smaller than 
that expected from the jet scenario. Furthermore, the mid-infrared spectrum obtained with the Spitzer Space Telescope indicates that 
the 10 jim silicate emission is optically thin and the silicate grains are of sub-micron size. This lends support to the presence of a 
circum-companion disk, because an optically thin circumbinary disk consisting of such small grains would be blown away by the 
intense radiation pressure of the primary (carbon-rich) star. If we assume Keplerian rotation for the circum-companion disk, the mass 
of the companion is estimated to be 0.5-0.8 Mq. We also identify CO2 emission features at 13-16 fim in the Spitzer spectrum of 
EU And — the first unambiguous detection of CO2 in silicate carbon stars. 
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1. Introduction 

Oxygen-rich circumstellar material (i.e., silicate and/or AI2O3) 
is usually associated with M-type asymptotic giant branch 
(AGB) stars, reflecting their photospheric chemical composi- 
tion. Surprisingly, however, carbon stars showing silicate emis- 
sion (so-called "silicate carbon stars") were discovered by IRAS 
(Little-Marenin I 19861 Willems & de Jong 1986). The subsequent 
detection of H2O and OH masers toward some of the silicate 
carbon stars, if not all, confirmed the oxygen-rich nature of the 
circumstellar material (e.g., Nakada et al. 119871 119881 Benson 
& Little- Mareninfl 987 ; Litde-Marenin et al. |1988l Bambaum et 
al. |19911 Engels 1994). On the other hand, optical spectroscopic 
studies show that '^C/'-'C ratios in silicate carbon stars are as 
low as 4-5 and thus they are classified as "J-type" carbon stars 
(e.g., Ohnaka & Tsuji 1999), which is diflicult to explain by the 
standard stellar evolution theory. 

One cuiTently accepted hypothesis suggests that silicate car- 
bon stars have a low-luminosity companion (possibly a main- 
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sequence star or a white dwarf) and that oxygen-rich material 
was shed by mass loss when the primary star was an oxygen- 
rich giant. The oxygen-rich material is stored in a circumbi- 
nary disk even after the primary star becomes a carbon star 
(Morris [1987t Llovd-Evans |199()b - Alternatively, Yamamura et 
al. ( 2000) propose that the oxygen-rich material is stored in a 
circumstellar disk around the companion. In this scenario, the 
observed silicate emission cannot originate directly from the 
circum-companion disk heated by the companion itself, because 
its luminosity is too low. Instead, Yamamura et al. (2000) argue 
that radiation pressure from the primary, carbon-rich, AGB star 
drives and heats an outflow from the circum-companion disk. 
The observed silicate emission originates from this outflow. 

Our recent A^-band (8-13 //m) spectro-interferometric obser- 
vations of the silicate carbon star IRAS08002-3803 using the 
ESO's Very Large Telescope Interferometer (VLTI) have spa- 
tially resolved the dusty environment of a silicate carbon star for 
the first time, and our radiative transfer modeling shows that the 
A^-band visibilities can be fairly explained by an optically thick 
circumbinary disk in which small (~0. 1 yum), amorphous silicate 
and a second grain species — amorphous carbon, large silicate 
grains (~5 jum), or metallic iron — coexist (Ohnaka et al. 2006J. 
The A^-band visibilities observed on another silicate carbon star 
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IRAS 18006-32 13 show a wavelength dependence very similar 
to 1RAS08002-3803, suggesting that IRAS 18006-32 13 also has 
a circumbinary disk (Deroo et al. 2007 ). 

Some evidence of long-lived disks around silicate carbon 
stars has been discovered by radio observations as well. Jura 
& Kahane d 19991 1 detected narrow CO (7=2-1, 1-0) emission 
lines toward two silicate carbon stars EU And and BM Gem, 
which indicate the presence of a reservoir of orbiting gas. 
Szczerba et al. ( I2006I I and Engels (priv. comm.) obtained high- 
resolution 22 GHz H2O maser maps toward the silicate carbon 
star V778 Cyg using MERLIN and VLBA, respectively. The 
H2O maser distributions are linearly aligned with a slightly S- 
shaped structure, which Szczerba et al. (2006) and Babkovskaia 
et al. ( I2OO6I 1 interpret as a warped circum-companion disk 
viewed almost edge-on. Therefore, as discussed in Yamamura et 
al. (2000) and Ohnaka et al. (2006), there may be two classes of 
silicate carbon stars: systems with optically thick circumbinary 
disks and those with circum-companion disks. 

While direct detection of companions of silicate carbon stars 
is still difficult due to the huge luminosity contrast (~ 10"* Lq 
and ^1 Lq for the primary star and the companion, respec- 
tively), Izumiura ( 120031 ) and Izumiura et al. ( 120071 1 detected 
blue continuum emission at ^4000 A and Balmer lines showing 
P Cygni profiles with an outflow velocity of ^400 km s ' toward 
BM Gem, which strongly suggests the presence of an accretion 
disk around an unseen companion. 

Despite this observational progress, the formation mecha- 
nisms of the circumbinary or circum-companion disks as well 
as the peculiar photospheric chemical composition of silicate 
carbon stars are little understood. Radio interferometry obser- 
vations of the maser emission produced by oxygen-bearing 
molecules provide us with an excellent opportunity to study the 
spatio-kinematic distribution of oxygen-rich gas around silicate 
carbon stars. In this paper, we present the results of multi-epoch, 
high-resolution observations of the 22 GHz H2O masers toward 
the silicate carbon star EU And. 



2. Observations 

The H2O maser emission associated with EU And 
(q-=23''19'"58.8814', 5=47°14'34.567", J2000.0, Epoch 2000.0, 
NOMAD Catalog, Zacharias et al. 2003F1 was observed at a 
rest frequency of 22.23508 GHz, using the Very Long Baseline 
Array (VLBA) operated by the National Radio Astronomy 
Observatory (NRAOfl The observations were carried out 
over three epochs spanning approximately 14 months: 2005 
November 5, 2006 March 19, and 2007 January 13. VLBAs 
10 antennas were used at the first and second epochs, while 
only 9 antennas could be used at the third epoch due to a 
technical problem of the St. Croix station. For each 5-hr epoch, 
EU And and two extragalactic calibrator sources (3C454.3 
and J2322-I-5057) were observed. We used the technique of 
rapid switching between the target source and the nearby phase 
reference source, J2322-I-5057, in order to remove residual 
phase offsets in the target source due to the atmosphere. The 



' This is offset by 32" from the position used for the VLA observa- 
tion on 1990 June 2 by Colomer et al. (2000), who detected no H2O 
maser toward EU And. The negative detection may have been due to 
this positional offset, although it is still possible that the H2O masers 
were indeed absent in June 1990. 

^ The National Radio Astronomy Observatory is a facility of the 
National Science Foundation operated under cooperative agreement by 
Associated Universities, Inc. 



data were recorded in dual circular polarization with the 8-MHz 
band centered on a local standard of rest (LSR) velocity of 
-36.0 kms"'. The system temperatures and the sensitivities 
were of the order of 100 K and 10 Jy K respectively, for all 
three epochs. 

The data were correlated at the VLBA coiTelator in Socoito, 
New Mexico, which produced auto- and cross-correlation spec- 
tra with 512 channels, corresponding to a channel spacing of 
15.63 kHz (0.22 kms '). For the calibration of the coiTe- 
lated data, we followed the standard reduction procedures for 
VLBA spectral line experiments using the Astronomical Image 
Processing System (AIPS) maintained by NRAO. A bandpass 
calibration was caiTied out using intermittent (every 30 minutes) 
scans on 3C453.3. A fringe fit was performed on the phase- 
reference source in order to remove effects due to instrumen- 
tation and atmosphere on the phase not removed in the correla- 
tion. The resulting residual phase delays, phase rates, and phases 
were applied to the target source. At this point a preliminary im- 
age cube was produced for EU And. Because the signal-to-noise 
in these images was not optimal, we performed a second fringe 
fit on a strong reference maser feature (at Vlsr - -34 km s"' ) 
in the spectrum of EU And and applied the residual phase rates 
to the target. Finally, an iterative self-calibration and imaging 
procedure was performed on the reference channel and the so- 
lutions for the residual phase and amplitude coiTections were 
applied to all spectral channels. Images with 1024x1024 pixels 
(61x61 mas) were produced for all velocity channels between 
Vlsr - -24 and -52 km s"' using beam sizes of 0.71x0.29 mas, 
0.73x0.36 mas, and 0.82x0.31 mas for the first, second, and 
third epochs, respectively. Typical rms off-source noise in the 
final images is 7-10 mJy beam"'. 

Maser components were identified, and their positions were 
measured by fitting the emission in the final images of each 
spectral channel with two-dimensional (2-D) Gaussians using 
the AIPS task SAD. The errors of the positions (within an 
epoch) are 10-20 //as for the strong maser features and approach 
~50 /ias for the weaker features. The absolute astrometric 
information obtained by phase-referencing is lost in the above 
self-calibration procedure. To restore this information to the 
final component positions, we fit 2-D Gaussians to the strongest 
maser feature (assumed to be the same for all three epochs at 
Vlsr = -34 km s"') prior to any self-calibration. The resulting 
shifts were applied to all maser components identified in the final 
images. The absolute positions of the strongest maser obtained 
for the first, second, and third epochs are (23''19'"58.8822\ 
47°14'34.550"), (23''19"'58.8823^ 47°14'34.548"), and 
(23''19'"58.8824', 47°14'34.546"), respectively (all in J2000.0). 

The accuracy of the absolute position of the strongest maser 
feature is primarily limited by errors resulting from the transfer 
of the residual phase corrections from the extragalactic refer- 
ence source to the target source and by the error in the absolute 
position of the reference source itself. In order to estimate the er- 
ror involved in the transfer of the residual phase coiTections, we 
essentially reversed the procedure used to determine the abso- 
lute position of the target maser source. We applied the residual 
phase solutions resulting from the self-calibration on the target 
source (strongest maser at Vlsr = -34 km s"' ) to the extra- 
galactic reference source. We then imaged the reference source 
and performed a 2-D Gaussian fit to determine its position. The 
separation between the positions of the target/calibratorpair (i.e., 
the arc length) was computed for 1) both sources calibrated with 
the extragalactic reference source solutions, and 2) both sources 
calibrated with the target source solutions. The difference be- 
tween the arc lengths yields an estimate of the error involved in 
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Fig. 1. 22 GHz H2O masers observed toward EU And with VLBA at three epochs. The positions and velocities of the H2O maser 
spots are derived by two-dimensional Gaussian fitting as described in Sect. |2] Left: Spatial distribution of the H2O masers with 
point-color representing the corresponding velocity bin in the spectrum and point-size proportional to the logarithm of the maser 
flux density. The origin of R.A. and DEC. offsets is (23''19"'58.8822', 47°14'34.550", J2000.0). The maser spots observed at three 
epochs are positionally offset because of proper motion and parallax. Right: Spectra formed by plotting maser component flux 
density versus LSR velocity, color-coded according to maser velocity. The solid lines represent vector-averaged cross-power spectra 
on the Los Alamos-Pie Town VLBA baseline. A color version of this figure is available in the electronic edition. 



the phase transfer The computation of the arc lengths is unnec- 
essary in experiments specifically designed for astrometry where 
a second extragalactic check source is typically observed. This 
procedure was repeated for all three epochs and the standard de- 
viation of the difference in arc lengths was found to be 0. 12 mas. 
This error in the phase transfer compares favorably with the error 
in the absolute position of the reference source itself, which is a 
defining source of the International Celestial Reference Frame 
(ICRF) with an eiTor of 0.26 mas (Fey et al. 2004). The total 
error in the absolute position of the strongest maser feature is 
estimated from the root-sum-square of the two errors described 
above to be ~0.3 mas. 

3. Results 

Figure [T] shows the spatio-kinematic distributions (left panel) 
and spectra (right panels) of the H2O masers toward EU And 
as observed over the three epochs. At all epochs, the H2O 
maser spectra are characterized by two strong peaks at Vlsr 
= -34 km s ' and -42 km s In the first and second epochs, 
weak features were also detected at -38 km s"', approximately 
in the middle of the two peaks. The velocity of these weak 
maser features is in agreement with that of the narrow CO lines 
( 37 km s"') measured by Jura & Kahane ("1999). The ve- 
locity of the weak maser features also coincides with the radial 
velocities derived from the optical spectra, which represent the 
velocity of the primary carbon-rich AGB star (Barnbaum |1991| ). 
However, as discussed below, the putative companion, instead 
of the primary star, is more likely to be located at the position 
of these weak maser components. We note that while the spectra 
taken over 14 months appear rather stable, the H2O masers to- 
ward EU And exhibit remarkable temporal variations in strength 



and velocity. For example, Little-Marenin ( 119881 ) detected maser 
components near Vlsr - -30 km s~^ as strong as 8 Jy, which are 
non-existent in our VLBA data. 

The H2O masers at all three epochs are linearly aligned 
across ~20 mas, with a slight hint of an S-shaped structure. 
As explained in Sect. |2] the offsets among the epochs shown 
in Fig. [T] are real and represent the proper motion and parallax. 
In fact, the slight curvature seen in the direction of the posi- 
tional displacement over the three epochs, which is discernible 
particularly for the strong features (the one at (0, 0) at the first 
epoch), shows the annual parallax component. The observed S- 
shaped structure can be interpreted as an edge-on disk or an 
emerging helical jet. Such a nearly linear distribution of the H2O 
masers is similar to those observed toward the so-called wa- 
ter fountain sources (e.g., Imai et al. 120021 Boboltz & Marvel 
120071 and references therein), which show well-collimated, fast 
jets of H2O masers. However, the outflow velocities of the H2O 
masers toward the water fountain sources are much higher (~60- 
150 km s"') than that observed toward EU And (~5 km s"'), al- 
though the projection effect is uncertain. 

In order to examine the jet scenario, we measured the in- 
crease of the separation (arc length) between the masers asso- 
ciated with the two highest peaks of emission in the spectrum. 
Since a single maser spans multiple channels in the spectrum, 
we performed a flux density squared weighted average over ve- 
locity and position prior to measuring the separations between 
the two peak features. The difference in the arc lengths between 
the first and third epochs is ~0.3 mas. Using this measurement, 
we estimate the velocity V and the inclination angle / (angle 
toward us out of the plane of the sky) of a water-fountain-like 
outflow for EU And as follows. The observed velocity separa- 
tion of the two strongest masers is 8 km s"', which means that 
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Fig. 2. The optical position of EU And in 2000 (NOMAD cat- 
alog) and the expected optical position at the first epoch of the 
VLBA observations (Nov. 2005). The error bars are the sum of 
the uncertainties of the absolute position and the proper motion 
(see Sect. Ell. 

2ysin! = 8 km s~'. The spatial separation of the two masers 
increases by 2y cos / x At/d = 0.3 mas, where At is the time in- 
terval between the first and third epochs (14 months), and d is the 
distance of EU And. With distances of 1 .5-2.6 kpc adopted (Jura 
& Kahane [19991 Engels [T994l l, we obtain V = 4.1^.3 km s ' 
and i = 68-77°. These outflow velocities are much lower than 
those observed toward the water-fountain sources, making the 
helical jet interpretation unfavorable, although the distance of 
EU And may be remarkably larger than the above estimates 
and/or the outflow velocity of EU And may indeed be so low 
due to some driving mechanism different from that in the water- 
fountain sources. A further long-term monitoring as well as a 
measurement of proper motion with VLBA would be necessary 
to settle this issue. 

The H2O maser distributions observed toward EU And are 
similar to that observed toward V778 Cyg by Szczerba et al. 
( 1200 6') and Engels (priv. comm.). Our observations of another 
sihcate carbon star IRAS07204-1032 using VLBA and the Very 
Large Array (VLA) have also revealed a similar spectrum and 
spatial distribution of the H2O masers (Boboltz et al. in prep.). 
The expectation of observing linear distributions of masers to- 
ward the three silicate carbon stars imaged to date may seem 
statistically unlikely, given a random distribution of disk incli- 
nation angles. However, it is quite possible that this is due to a 
selection effect in which masers are preferentially detected in ob- 
jects with the most favorable geometry for maser amplification 
(i.e., nearly edge-on disks). A survey of more objects with VLA 
and VLBA would be useful for obtaining a definitive answer to 
the interpretation of the masers toward silicate carbon stars. 

4. Discussion 

Although our H2O maser maps of EU And lend support to 
the presence of an edge-on disk, these observations alone can- 
not clarify whether the disk is located around the carbon- 
rich primary star (circum-primary disk) or the whole binary 
system (circumbinary disk) or an unseen companion (circum- 
companion disk). The optical position of EU And, which corre- 



sponds to that of the primary stai-, is (a, 6) = (23''19'"58.8814^ 
+47°14'34.567", J2000.0) with proper motion of (//„, i^g) = 
(5. 5 + 1. 7, -1.5 + 2.0) mas yr"' (NOMAD catalog, Zacharias et 
al. 120031 1. The position of the primary star expected for the first 
epoch, which is shown in Fig. |2l appears to be offset from the 
H2O masers. However, the uncertainty of the absolute position 
(+15 mas, Zacharias et al.'2004) as well as that of the propermo- 
tion is too large to examine whether the primary star is located 
at the center of the H2O masers or indeed offset from the maser 
distribution. 

4. 1 . Constraints from the mid-infrared spectrum 

On the other hand, as Yamamura et al. (I2000I I argue, a circum- 
primary or circumbinary disk would be unstable against the 
intense radiation pressure from the primary carbon-rich star 
(~10^-10'* L0), if the disk is entirely optically thin. In order 
to examine whether or not the silicate dust around EU And 
is optically thin, we obtained its mid-infrared spectrum from 
the data archive of the Spitzer Space Telescope (Werner et al. 
,2004) . EU And was observed on 2004 December 9 (Program 
ID: P03235, P.I.: C. Waelk ens) with the InfraRed Spectrograplfl 
(IRS, Houck et al. '2004| in the Short-High (SH) and Long- 
High (LH) modes with a spectral resolution of ~600 . We down- 
loaded the Basic Calibrated Data (BCD) processed with the S 15- 
3 pipeline and extracted the spectrum using SMART v.6.2.5 
(Higdon et al. |2004l Pi Figure [3l shows the Spitzer/IRS spec- 
trum of EU And (no sky subtraction was performed) together 
with the spectra of the silicate carbon star V778 Cyg and the 
oxygen-rich AGB star o Cet obtained with the Short Wavelength 
Spectrometer (SWS) onboard the Infrared Space Observatory 
(ISO). As Yamamura et al. (120001 ) show, the 10 /im silicate 
features of the latter two stars indicate that the silicate emis- 
sion is optically thin and silicate grains are of sub-micron size. 
Obviously, the spectrum of EU And closely resembles those of 
V778 Cyg and o Cet. A comparison with the IRAS 12 jjm flux of 
EU And and its Spitzer/IRS spectrum reveals that the flux level 
has little changed for the last 21 years, and EU And is similar 
to V778 Cyg in this aspect as well. This means that the sili- 
cate emission of EU And is optically thin with sub-micron grain 
sizes, rendering the possibility of a circum-primary or circumbi- 
nary disk unlikely. Therefore, the disk-like structure discovered 
by the H2O masers is likely to represent a circum-companion 
disk. 

The predominance of sub-micron-sized grains appears to 
contradict the conclusion of Jura & Kahane (1999) that the grain 
size around EU And should be as large as >0.2 cm to be grav- 
itationally bound against radiation pressure. However, they as- 
sumed a circum-primary or circumbinary disk around the cen- 
tral star with a luminosity of 10"* Lq. In the case of a disk 
around a low-luminosity companion (~1 Lq), grains as small as 
0.2 yum are stable against the radiation pressure of the compan- 
ion, although this does not entirely exclude the presence of larger 
grains due to grain growth in dense regions (e.g., mid-plane) of 
the circum-companion disk. 

In the Spitzer/IRS spectrum of EU And shown in Fig. |3l we 
can identify the CO2 emission features. To illustrate this identi- 
fication, we also show the (scaled) absorption cross sections of 



^ The IRS was a collaborative venture between Cornell University 
and Ball Aerospace Corporation funded by NASA through the Jet 
Propulsion Laboratory and Ames Research Center. 

SMART was developed by the IRS Team at Cornell University and 
is available through the Spitzer Science Center at Caltech. 
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CO2 and CO2, which were calculated using the CO2 line list 
of the HITRAN database (Rothman et al.'2005^ for 1000 K (the 
column densities of the both isotopic species are set to be equal). 
The CO2 features are clearly seen in the ISO/SWS spectrum of 
the oxygen-rich star o Cet, and also marginally in V778 Cyg. 
Yamamura et al. (120001) could only tentatively identify these CO2 
features in V778 Cyg because of the low S/N ratio, and the case 
of EU And is the first unambiguous detection of CO2 toward sil- 
icate carbon stars. This is another piece of evidence for oxygen- 
rich gas around silicate carbon stars, and the mid-infrared CO2 
features will enable us to estimate the CO2 density and tempera- 
ture if the disk structure is more tightly constrained in the future. 
It would also be possible to derive the '^C/'^C ratio of the cir- 
cumstellar gas from an analysis of the mid-infrared ^^C02 and 
'^C02 features, which would shed new light on the origin of the 
abnormally low '^C/'^'C ratios in silicate carbon stars. 

It is worth noting that o Cet (Mira AB), whose mid-infrared 
spectrum closely resembles that of EU And, is a well-known bi- 
nary system (separation ~0'.'5) with a main-sequence or white 
dwarf companion. The high-resolution mid-infrared images of 
Mira AB recently obtained by Ireland et al. ( |2007| l suggest that 
the edge of an accretion disk around Mira B (companion) is 
heated by Mka A (the primary star at the AGB). As Izumiura 
et al. ( 120071 ) suggest, this may resemble the circumstellar envi- 
ronment of some — if not all — silicate carbon stars (including 
EU And). Unlike silicate carbon stars, the circum-companion 
disk around Mira B consisting of silicate is embedded in the 
oxygen-rich outflow from Mira A, which masks the mid-infrared 
spectroscopic signature of the circum-companion disk. 

4.2. Nature of the companion 

In order to draw a p-V diagram for the observed masers to- 
ward EU And, we determined the center of the overall spa- 
tial distribution by taking the mid-point of a line connecting 
the mean positions of the redshifted and blueshifted masers at 
each epoch. The center velocity was also derived by halving 
the difference between the mean velocities of the redshifted and 
blueshifted masers. The resulting p-V diagram for EU And is 
shown in Fig. |4] It is very similar to that derived for V778 Cyg 
by Szczerba et al. (120061) . although the number of maser spots 
is much smaller than for V778 Cyg. Both objects show a lin- 
ear part with very large gradients at both ends of the veloc- 
ity range. This is consistent with the presence of a Keplerian 
disk, although it is not definitive evidence. If we assume that the 
(circum-companion) disk is in Keplerian rotation, and the center 
of the disk is located in the middle of the maser distribution at 
Vlsr = -37.9 km s"', we can estimate the mass of the compan- 
ion at the disk center. As shown in Pestalozzi et al. (2004), maser 
emission observed toward an edge-on disk has three maxima — 
on the line of sight toward the disk center and near the lines of 
sight tangential to the outer edge of the disk. If the strong H2O 
maser components to NE and S W correspond to the maxima near 
the outer edge of the disk, it follows that the Keplerian velocity 
is ~5 km s"' at a radius of 10 mas = 10 x (d/kpc) AU, where d is 
the distance of EU And. This translates into 0.3 (d/kpc) Mq for 
the central mass (0.5-0.8 Mq for d - 1.5-2.6 kpc). This mass 
would be reasonable for an unevolved main-sequence star or a 
white dwarf. 

It should be noted, however, that the strong maser compo- 
nents at -42 and -34 km s"' may not represent the outer edge 
of the disk. As mentioned above, the strong maser components at 
~ -30 km s"' were detected by Litde-Marenin et al. ( 119881 ), but 
the lack of spatial information does not allow us to take these ob- 
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Fig. 3. a: Spitzer/IRS spectrum of EU And, together with the 
scaled ISO/SWS spectra of V778 Cyg and o Cet. The diamond 
represents the IRAS 12 fim flux of EU And. b: CO2 emission 
features identified in the Spitzer/IRS spectrum of EU And are 
marked with the ticks. The scaled absorption cross section of 
'2C02 and ^^C02 for 1000 K are also shown. 

servations into account in estimating the mass of the secondary 
star. Moreover, there may be masers fainter than the VLBA's 
detection limit and/or more diffuse emission which would be re- 
solved out by VLBA. Therefore, the above estimate of the com- 
panion mass based on Keplerian rotation is a lower limit. The 
estimate of the companion mass also depends on the disk geom- 
etry assumed in the analysis. For example, for the MERLIN data 
on V778 Cyg presented by Szczerba et al. (I2006I I. Babkovskaia 
et al. ("2006) applied a doubly warped disk model and derived a 
companion mass of 1.7 Mq, which is significantly larger than 
the 0.06 Mq derived by Szczerba et al. (l2006l l. This highlights, 
once again, the importance of long-term monitoring of the H2O 
masers in order to derive the properties of the disk and the sec- 
ondary star 

5. Conclusion 

Our VLBA observations of the 22 GHz H2O masers toward 
the silicate carbon star EU And have revealed that the masers 
are linearly aligned with a slight hint of an S-shaped structure, 
with the masers in SW redshifted and those in NE blueshifted. 
Such a spatio-kinematic structure can be interpreted either as 
an edge-on disk or an emerging jet. The H2O maser maps ob- 
tained at three epochs over 14 months show a little outward 
motion of 0.3 mas, but the outflow velocity of ~4 km s"' esti- 
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Fig. 4. The p-V diagram for the H2O masers toward EU And. 
The center position and velocity were derived as described in 
Sect. 14.21 The center position (for the first epoch) is located at 
(7.7, 5.9) mas in the left panel of Fig. [1] while the center velocity 
Vc is -37.9 km s"'. The circles, squares, and triangles represent 
the data taken at the first, second, and third epochs, respectively, 
with the velocities represented by the colors as in Fig.[T] A color 
version of this figure is available in the electronic edition. 
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mated from this outward motion is too low compared to those 
observed toward the water-fountain sources. This lends support 
to the disk interpretation, although some kind of outflow (pos- 
sibly with a velocity much lower than in the water-fountain 
sources) cannot be entirely ruled out due to the uncertainty in 
the distance. The mid-infrared spectrum of EU And obtained 
with the Spitzer/IRS shows that the silicate emission is optically 
thin and emanates from sub-micron-sized grains, which suggests 
that the H2O masers originate in a circum-companion disk seen 
nearly edge-on. Furthermore, we unambiguously identified the 
CO2 features at 13-16 //m in the Spitzer spectrum for the first 
time. 

If we assume that the disk is in Keplerian motion, the mass 
of the putative secondary star is estimated to be 0.5-0.8 Mq. 
However, given the remarkable variability of the water masers 
and the possible presence of fainter or more diffuse emission not 
detected by VLBA, the estimated mass should be regarded as a 
lower limit. Observations with a more compact, high-sensitivity, 
array (e.g., the VLA or MERLIN) might be useful in character- 
izing any such extended emission. Future VLBA observations 
will be necessary to put stronger constraints on the geometry of 
the disk and the properties of the secondary star. Astrometric 
monitoring with the VLBA should also provide more accurate 
estimates of the parallax and proper motion of EU And, thus im- 
proving our understanding of the three-dimensional structure of 
the oxygen-rich gas surrounding the system. 
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